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Abstract. Dense star clusters expand until their sizes are limited by the tidal field 
of their host galaxy. During this expansion phase the member stars evolve and lose 
mass. We show that for clusters with short initial relaxation time scales (<100 Myr) 
the dynamical expansion is largely powered by mass loss from stars in the core, but 
happens on a relaxation time scale. That is, the energy release following stellar mass 
loss is in balance with the amount of energy that is transported outward by two-body 
relaxation. 



1. Introduction 

Like a star, a stellar cluster produces heat by contracting. The positive energy that is 
released by the more tightly bo und core is diffused outward by two- body relaxation 



dLvnden-Bell & Eggleton|[l980l) on a relaxation time scale (equation [T] Spitz er & Hartl 



19711) . Just as for a star, th e cont raction phase in a cluster is only a small fraction 



of the total lifetime^. HenorJ ( 1975b had the insight that the evolution of clusters after 



core contraction can be understood by assuming that an energy source provides the right 
amount of energy in a self-regulating fashion to sustain the two-body relaxation process. 
A breakthrough in the understanding of cluster evolution, comparable to Eddington's 

work on stellar structure. 

In single-mass clusters the energy is s upplied by binary star s dHenonlll96lb that 



form in the densest phase of core collapse dTanikawa et al.ll2012f) . They provide en- 



ergy by int eracting with other stars, thereby increasing the total (negative) energy of 
the cluster (Heggid ll975h . In more realistic clusters mass loss of member stars also in- 
creases the total (negative) energy. In this contribution we show that the corresponding 
energy production rate can come into balance with the flow of energy due to two-body 
relaxation. 



2. Time scales 

When considering mass loss of stars in a collisional system there are a few time scales 
that need to be considered: the timescale for massive stars to segregate to the cluster 
core by dynamical friction, and the lifetime of stars, or main sequence time, r ms . The 
former depends on the mass of the star, m, the mean stellar mass, (m) and the half-mass 



'This is true for clusters that are initially much denser than the density within the tidal radius. 
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relaxation time scale of the cluster, r r h, as tm - (m/(m) ) 1 t± (Sp itzer|[l969h . We use 



the conventional expression for r r h ( Spitzer & H artHl97l1) 



N V2 3/2 

r rh = 0.138—^2 . (1) 

^/G(m)lnA 

Here ./V is the number of stars, is the half-mass radius, G is the gravitational con- 
stant and In A is the Coulomb logarith m. The argument of th e Coulomb logarithm 
is A - 0.1 IN for single-mass syst ems dGiersz & Heggiel fl994h and A - 0.02N for 
clusters with a wide mass spectrum dGiersz & Heggielll996l) . The mass-loss rate is pre- 
dominantly set by r ms and this timescale depends on the properties of the star, mainly 
m, and not on the dynamical state of the cluster. The two timescales r^f and r ms be- 
come comparable for the most massive main sequence star after some period of time. 
This is because r ms « m~ A , with 2.5 < A < 3, and Tdf °c m~ l . A star for which 
Tdf(m) = r ms (m) reaches the core within its lifetime and will lose its mass there. From 
that moment onwards mass loss will proceed on average from the central parts of the 
cluster. The amount of energy that is released depends on the mass and specific energy 
of the star at the moment it loses mass. In § [3] we quantify how much energy a cluster 
needs to evolve and in §@]we show how much energy can be generated by mass loss. 



3. Energy demand 

Let us first consider how much energy is needed per unit of time. Just as for a star, 
the flow of energy, or luminosity E, is determined by the properties of the syste m as a 



whole. It equals a fraction C of the instantane ous total energy E per unit of T r h (IHenon 
119611 : iHeggie & Hull2003l : iGieles et al.ll20ll 

E C 

-5 = -—■ (2) 
The total energy for a cluster in virial equilibrium may be written as 

GM 2 

E = -k , (3) 

Hi 

where k - 0.20 - 0.25, depending on the density profile and M = (m)N is the total 
mass. The evolution of the main cluster parameters (Af and r\) can be understood once 
we know how r r h evolves. For this it is insightful to consider the idealised models of 
Henon. In the absence of a tidal field the energy flow leads to a self-similar exp ansion 
at a r oughly constant mass and r r h grows linearly in time such that x t 2 ^ (|Henonl 
Il965h . In a stea dy tidal field the radius is limited by the tidal radius and r r h decreases 



linearly in time (lHenor][l961h . In b oth cases £^0.1 (see the d erivation from Henon's 



work in Giel es et al.l 1201 ll and see lAlexander & Gielei |2012| for measurements of £ 



from numerical simulations) and £ can be more than a n order of magnitude higher for 
clusters with a wide mass spectrum dGieles et al.ll2010h . 

We consider clusters that are not strongly confined by a tidal field. This is justi- 
fied by the observation that about half of the (old) Galactic g l obular clusters are still 
expanding towards their tidal boundary dBaumgardt et al.ll2010l : lGieles et al. 2011). We 
then find from equation (Q]) and Q that r r h = (3/2)£? and 
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This means that the (absolute value of the) cluster energy decreases by roughly 80% 
each age dex during the expansion. It is in this phase that stellar mass loss can provide 
the energy. Although a cluster only loses about half of its initial mass due to stellar 
mass loss over a Hubble time, this modest mass-loss rate in fact generates energy at a 
rate comparable to what is required (equation |4]), as we will show in §|4] 



4. Energy supplied by mass loss 

As we aim to show in this section, the rate of energy increase as a result of mass loss 
of stars, E&, leads to a power-law decline of —E (equation |U). This is primarily because 
the evolution of the total mass as a result of stellar evolution from a stellar system 
with a Salpeter like stellar initial mass function, declines roughly as a power of time 
M(f) oc r y , with v 0.07 (this corresponds to a loss of about 15% of the mass every 
age dex, e.g. iBruzual & Charlot]|2003l) . The instantaneous mass-loss rate is thus 

^ = - V -. (5) 
M t 

To relate this to energy we need to know from where and how fast mass is lost. 
We assume that the mass-loss time scale is much longer than the local c rossing time , 
such that the cluster responds adiabatically and retains virial equilibrium. iHills! ( 1980h 



discussed the adiabatic response of a gravitational system if mass is lost homologous 
with the density profile of the cluster, that is, without a preferred location. The cluster 
radius then increases as M" 1 and the change in energy relates to the change in mass as 
(equation [3]) 

dE dM 

T - ~M' (6) 

When Tdf << r ms mass loss happens predominantly from the core of the cluster, 
because the most massive stars, i.e. the ones contributing most to the total stellar mass 
loss, are centrally concentrated because of dynamical friction^- Mass loss affects both 
the potential energy U = 2E and the kinetic energy T = (l/2)M(v 2 ), where (v 2 ) is the 
mean-square velocity of the system. This is because the mass that is removed carries 
away a bit of potential and kinetic energy and because of the subsequent redistribution 
of potential and kinetic energy of the remaining stars. We can express the change in the 
potential energy in terms of the specific potential of the mass that is lost, 0a . relative to 
the average specific potential of all the stars, (<p) = 2U/M, as dU/U = 2(0 A /(0»dM/M. 
For the kinetic term we need to compare the squared velocity of the mass that is lost, 
v 2 A , to the mean square velocity of the cluster: dT/T = {v 2 J (v 2 ))dM I M . Using our 
assumption that the remaining stars quickly restore virial equilibrium (E = U/2 = -T) 
we can relate the change in the total energy to the change in mass 



2 This simple view can get slightly more complicated if the turn-off mass is significantly lower than the 
non-evolving remnants. If many remnants are retained they will eventually dominate the dynamics in the 
inner regions of cluster. 
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** = 2^-* (7) 
E U T w 

= x ™, with 4^-4- (8) 

Af <0> (v 2 > 

For homologous mass loss (0a) = (0) and (v 2 ) = (v 2 ) and therefore x = 3. This is 
Hills' result (equation [6]). The general result for the rate of energy 'supply' is 

^ = --• (9) 
E t 

From a comparison to the energy 'demand' (equation |U) we see that a balance 
between energy supply and demand is reached if xv - 2/3. Because v — 0.07 the con- 
dition for energy generation by stellar mass loss to be in balance by two-body relaxation 
becomes x - 10. To answer the question whether this value is realistic, we need to con- 
sider a cluster poten t ial tha t described a cluster in the balanced evolution phase. One 



example is the JaffeJ (119831) model, which has the required r density cusp. From the 
potential of this model we find <p(r)/(<f>) = ln(l + r^/r). If we ignore the small contri- 
bution of the kinetic term to x in equation © we find that x - 10 if mass loss takes 
place at r - O.lrh. This is probably realistic in the case of a mass segregated system. 
This shows that the product xv can be of sufficient magnitude (about 2/3, equation |4]) 
to drive the relaxation process. In a forthcoming paper (Gieles, Heggie & Church) we 
will illustrate this result with numerical simulations and show that stellar mass loss is 
indeed a viable energy source in driving the dynamical expansion of star clusters. 
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